Preprint typeset in JHEP style - HYPER VERSION 



Hyperscaling- Violation on Probe D-Branes 



Martin Ammon", Matthias Kaminski , Andreas Karch 

° Department of Physics and Astronomy, University of California, Los Angeles, CA 
90095, USA 

■ , ^ Department of Physics, University of Washington, Seattle, WA 98195, USA 

(jj Email: ammon@physics.ucla.edu, mski@uw.edu, akarch@uw.edu 

Abstract: For the field theories dual to D3/D7- and D3/D5-brane systems we find non- 
relativistic finite density fixed points exhibiting a violation of hyperscaling. This violation 
is measured by the critical exponent 9 = 1 while the dynamical critical exponent is z = 2. 
At zero temperature we compute the thermodynamic potentials, the speed of normal sound, 
and the speed of zero sound for both these massive D3/D(2n+l)-brane systems near their 
non-relativistic fixed points. Moreover, we determine the first correction to the free energy 
for small temperatures yielding the critical exponents a and v. 
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1. Introduction 

Critical exponents traditionally arise in the context of dynamical critical phenomena ^ . They 
characterize the behavior of systems near a phase transition, organizing them into universality 
classes. These critical exponents satisfy a number of relations among each other. One of these 
is the hyperscaling relation ni' = 2 — a. The assumption underlying the hyperscaling relation 
is that all dimensionful quantities scale with their natural power of length. In particular, 
the entropy density should scale as a volume. In systems where the only length scale is 
the correlation length, this in particular implies that the entropy density should scale like 
the correlation length to the n, where n is the number of spatial dimensions. This relation 
can be violated in certain cases, and this violation is measured by the hyperscaling violation 
exponent 6 in the following way Q 

{n-e)i^ = 2-a, (1.1) 

where a describes the scaling of the specific heat capacity C ~ t~", while describes the 
scaling of the correlation length ^ ~ \t\~'^ near the phase transition at temperature T = Tc- 
n are the number of spatial dimensions and r is the reduced temperature given by r = 
(T — Tc)/Tc. Such theories show conformal symmetry at the transition T = Tc. However, the 
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scaling behavior in this context is always extracted from T ^ Tc, where the theory should be 
slightly deformed away from conformal invar iance. 

In mean field theory, critical exponents only depend on the form of the effective potential 
and so are independent of spatial dimension n. As such, mean field theory and hyperscaling 
are only consistent in a critical dimension njj, the so-called upper critical dimension. For the 
Ising model this critical dimension is njj = 4. The expectation is that for n < njj hyperscaling 
is valid, but mean field theory is not, whereas in n > njj mean field theory is valid and hence 
hyperscaling violated. As theories with a supergravity dual are intrinsically classical and 
often reproduce mean field results, one should presumably not be surprised that hyperscaling 
violation is actually the norm in these systems, as seen for example in @. 

So far we discussed thermal phase transitions occurring at Tc ^ 0. In the present paper, 
we consider a phase transition at zero temperature which is therefore driven by quantum 
fluctuations. Such quantum phase transitions occur at a particular value rc of the control 
parameter r, which can be a magnetic field, pressure - or in the case of the present paper a 
chemical potential fi. In the case of a quantum phase transition the relations stated above 
hold if we replace the number of spatial dimensions, n, by n + z. Here z is the dynamic 
critical exponent associated with the scaling of dynamic quantities. One such quantity is the 
correlation time tcorr oc which represents the characteristic time scale of the dynamical sys- 



tem under consideration. Now, for example, the hyperscaling relation (1.1) with hyperscaling 
violating exponent 6 reads 

{n + z-e)u = 2-a. (1.2) 

At non- vanishing temperature, T ^ 0, the free energy density / in the quantum critical region 
can be written as (see e.g. |^]) 



where 5 = r — rc- The function g is called scaling function and satisfies g{T = 0) = 1. 

Also for quantum phase transitions a critical dimension njj exists analogous to the clas- 
sical phase transitions mentioned above. But for quantum phase transitions we now have to 
consider the modified number of spatial dimensions. That is, we expect that for n + z < nu 
hyperscaling is valid, mean field theory is not. And vice versa in the case n + z > nij. 

Hyperscaling violation has been recently analyzed |^ ^ holographically In these 
works a scale was introduced by studying the theory at finite temperature ^; in the present 
paper we will similarly introduce a scale via a finite chemical potential. Neither in ||6|, |^ nor 
in the present paper does the static coefficient v mentioned above appear explicitly. However, 
we will see that in the present paper by construction vz = 1. The assumption of hyperscaling 
still should mean that all dimensionful quantities have their natural length dimension. This 
time the natural length scale is ~ T^l^ if we are studying the theory at finite temperature 
but zero density or ~ /i^/^ when the chemical potential is non-zero but the temperature 



^At small temperatures, the quantum critical region around r = rc is still governed by the critical point at 
= Tc and at Tc = 0. 
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vanishes. The entropy density in such a theory is given by s ~ r"/^ if hyperscahng is intact, or 
by s ~ fi"^-^)/^ if hyperscaling is violated P|. Roughly speaking, in theories with hyperscahng 
violation the length dimension of a spatial volume is reduced from n to n — 9. 

In general the gauge/gravity correspondence relates quantum field theories to gravity 
theories defined on spacetimes with metrics which are asymptotically Anti-deSitter (AdS). 
In Q a particular behavior of the metric on the gravity side was identified to be yielding 
non-trivial 9 and z on the gauge side. This particular behavior amounts to a specific scaling 
of the infrared (IR) limit of any AdS metric 

ds^ = p2 (^_p2n(z-l)/{n^e)^^2 ^ ^-2[2n-B) / (n-O) ^^2 ^ ^^.2^ ^ (^4) 

with the time direction t, n spatial directions Xj, and the emergent holographic direction 
p. From a rescaling with constant A we get t — )• A~^t, Xi — )• A^^Xj, p — )• A*^"~^)/"/9, and 
most prominently ds — )• A~^/"ds. At finite temperature this metric will develop a horizon at 
p = Ph- The horizon area is of size p^ and is proportional to the entropy of the dual field 
theory. As argued in P, we obtain a scaling of the entropy density s ~ — >• A"'~^s. This 
correspondence was further developed in [|, |, ||, . 

In this paper we consider a quantum field theory at zero temperature, which undergoes 
a quantum phase transition at a critical chemical potential pc- Both the D3/D7-, as well as 
D3/D5-brane systems have a dual description in terms of a gauge theory coupled to matter 
fields in the fundamental representation. In the quantum field theory these fundamental 
flavors have mass M, so at zero temperatures a chemical potential p < pc = M does not lead 
to any finite density, whereas for chemical potentials p > M a finite charge density develops. 
As shown in these two phases are separated by a second order phase transition. In the 
holographic dual, the transition corresponds to a topology change in the embedding of the 
probe branes (from a configuration where the probe branes are separated from the horizon 
to one where they cross the horizon). That is, in our system the chemical potential p acts 
as a tunable external parameter (analogous to the temperature in a ferromagnet). The order 
parameter is the charge density d (analogous to the magnetization density in a ferromagnet). 
The mass M of the fundamental fiavor particles provides a microscopic scale in this system 
determining the critical value of p (like the lattice constant and spin-couplings in the case of 
the ferromagnet). Despite the appearance of the microscopic scale M our system is described 
by a scale invariant quantum critical theory at the 2nd order phase transition point when 
p = pc- The non-relativistic chemical potential p = p — pc deforms the theory away from 
the critical point pc in the spirit of the statistical mechanics treatment reviewed above. Our 
main result is that the critical point at p = pc is governed by a scale invariant theory with 
non-trivial dynamical exponent z and a non-zero hyperscaling violation exponent 9, which 
acquire the values 

z = 2, 9 = 1. (1.5) 

This paper is structured as follows: We begin by sketching our line of argument in 
section ^. This leads to our main results, the values of the scaling exponents z and 9. In 
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section y we describe our holographic model, and prove details about its thermodynamics, 
normal and zero sound modes in the ensuing section ^ A generalization to finite temperature 
is carried out in section ^. We discuss our results in section |6|, while the appendix shows some 
details of the zero sound mode calculation. 



2. Scaling exponents and the main argument 

Consider a system with hyperscaling violation exponent and dynamical critical exponent z 
in n spatial dimensions. In such a system we work out scaling dimensions of thermodynamic 
quantities in all generality. 

Let us further consider a field theory at zero temperature, which only has one dimen- 
sionful scale, the chemical potential fl. The chemical potential p, can be written as the time 
component of a vector field and thus scales like an inverse time-coordinate with [p] = z. From 
thermodynamics at zero temperature we obtain the energy density e, and the grand canonical 
potential density is given hy 0, = e — fid in terms of the energy density e, charge density d. 
So since [E] = z, the energy density has [e] = n — 6 + z, because by definition the spatial 
volume scales like [V] = [d'^x] = — {n — 9). Thus, [e] = n — 9 + z and also [0] = [e] = n — 9 + z. 
Therefore the charge density has to have scaling dimension [d\ = n — 9. Up to a sign the 
pressure p is given by the density of the grand canonical potential and therefore we also know 
= n — 9 + z. At zero temperature the only dimensionful scale is /2 and we hence have to 
have 

p = -Vl = CoA("+"-^)/" (2.1) 

with some constant Cq. This implies that the energy density is 

n — 9 , ^ 

e = -p + fld= p, (2.2) 



where we have used d = dp/dfi. Note that (2.2) is a generalization of the commonly postulated 



ze = np \ n , 18] 



Heating up this theory to finite temperature T with [T] = z, we also obtain the scaling 
behavior of the entropy density [s] = n — 9. 

The goal of this paper is to realize a system as it is described above. Thus we consider 
a particular field theory: = 4 super- Yang-Mills theory in 3 + 1 dimensions with N = 2 
fundamental hypermultiplets, which all have the mass M. These hypermultiplets can live in 
either 3+1 or 2-|-l dimensions for the cases we consider. In addition our theory enjoys a U{1)- 
symmetry yielding a conserved charge, for which we switch on a density d, thermodynamically 
conjugate to the chemical potential ^. 

For fi = M both systems have a second order phase transition at zero temperature. In 
this paper we consider the regime /i ~ Af for /i > M. As an expansion parameter we use 
fl = fi — M <^ M, and consider only the leading non-vanishing contributions in fi/M. M is 
the rest mass of the particles in our hypermultiplets, and the chemical potential is the energy 
which is needed to add one such particle to the system. Therefore the limit p, <^ M can be 



-4- 



thought of as a non-relativistic Umit for this theory. The non-relativistic energy-density is 
given by 

e = e — ercst = —p + ^J-d — Md = —p + fid , (2-3) 

where e is the relativistic energy density, erest = Md the rest mass energy density of particles 
with mass M. 

In this paper we consider cases in which the hypermultiplet particles can propagate in 
3 + 1 and 2 + 1 dimensions (corresponding to D3/D7 and D3/D5 brane setups, respectively). 
Now, we state our first result about the non-relativistic energy density: 

P 

e = p , and ^ ~ 2 ' C^-^) 
for the 3 + 1- and the 2 + 1-dimensional case, respectively. These relations will be derived 



from D-brane thermodynamics below. Using (2.4) in ( |2.2D immediately yields the following 
relation between 6 and z: 

3-6 = z, and 2-6 = ^, (2.5) 

for 3 + 1 and 2 + 1 dimensions, respectively. 

Now, our second result states that the dispersion of the normal sound mode in both the 
3 + 1- and 2 + 1-dimensional theory is 

--/ffc, (2.6) 

for <C M. We will see later that this is confirmed also by the dispersion of the zero sound 
mode. The scaling dimension of the left-hand side is given by [u] = z, whereas the right-hand 
side gives [p]/2 + [k] = z/2 + 1. Therefore we conclude that z = 2. Using ( ^.5|) we obtain 
for both cases the value 6 = 1. Thus we find non-trivial dynamical scaling z and a non-zero 
hyperscaling violation as claimed in the introduction in (|1.5|). As we will see, mean field 
theory in this system yields z = 2 and 6 = n — 2. Hence, for the n = 3 case our findings 
are consistent with mean field theory. However, in n = 2 mean field theory would yield a 
vanishing hyperscaling violating exponent whereas we find again 9 = 1. 

Now let us prove case by case the two results we have used above to derive the scaling 
exponents in these cases. 

3. D3/D(2n+1) brane embeddings 

Let us now discuss the gravity dual of the field theory considered in section ^. The A/" = 4 
supersymmetric vector multiplet is realized by Nc coincident D3-branes which are aligned 
along the directions 0, 1,2,3 in the fiat ten dimensional spacetime. The Nf hypermultiplets 
propagating in n+1 dimensions (where n is the number of spatial dimensions) are represented 
by Nf D(2n+l)-branes which are aligned as shown in figure |l]. Concretely, we discuss 



the n = 3 D3/D7- |Tg and the n = 2 D3/D5-system [|g, |I| in what follows. 
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Figure 1: Coordinate directions in which the Dp-branes extend are marked by x. While D3- and 
D7-branes share the four Minkowski directions, D3- and D5-branes share only a 2 + 1-dimensional 
plane in the four Minkowski directions. 



Gravity setup Taking the large Nc limit AT^ — >■ oo, A = gyM^c fixed, as well as the strong 
coupling limit A S> 1 we can replace the Nc coincident D3-branes by their near-horizon 
geometry. In the limit Nf <^ Nc, the Nf D(2n-|-1) branes can be viewed as probe branes in 
this geometry. 

The background metric which we take to be AdS5 x can be written as 

ds'^ = H{r) r/^^ dxi'dx'' + H'^{r) (^p^ + p'^dsln + dy'^ + ^^dzf^ > (3-1) 
where H{r) = r'^/R'^ with 

4-n 

r' = p' + y' + ^zl (3.2) 

i=l 

and we have arranged the coordinates of in terms of an S'^ (n=2 for D3/D5 and n=3 
for D3/D7) as well as coordinates y and Zi (where i = 1,...,4 — n). These coordinates are 
very convenient to describe the embedding of the D(2n+l)-branes into S^. In particular, the 
D(2ra+l)-brane wraps the S"^, extends along the radial direction p, the time direction as well 
as 3 — n spatial directions. The world volume coordinates are denoted by ^. The embedding 
of the D(2n-|-l)-branes into AdS^ x is given by the transverse scalars y, Zi and x"' (where 
o, = n + 1, . . . , 3) as a function of the world volume coordinates ^. The embeddings at hand 
should neither break the rotational invariance of S*", nor the Lorentz-invariance of the common 
directions along 0,1,2,3. Therefore the embedding is specified by y = y{p), Zi = Zi{p), and 
= x"(p) where a = n + 1, . . . ,3. Here we will set Zi{p) = x"'{p) = 0. Besides the scalar 
fields, the dynamics of the probe D(2n+l)-branes is further described by a ?7(iVj;)-valued 
gauge field. Here we will restrict ourselves to the U{1) C U{Nf) gauge group and consider a 
non- vanishing gauge field component At{p). 

The dynamics of the D(2n+l)-brane is given by the DBI action 



SD{2n+l) = -NfTD^2n+l) J cfi""^^^^ -det{gab + It, a' Fab) , 



(3.3) 



where gab is the induced world-volume metric on the world volume of the D(2n+1) brane and 
Fab is the field strength tensor of the U(l) gauge field. Moreover 7£)(2n+i) is the tension of the 
D(2n+l)-brane. Since none of the fields depends on the coordinates of S"" or the field theory 
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directions in common, we can perform the integration over nine directions to get the action 
density, SD(2n+i) which is the action divided by the volume of the field theory directions in 
common 

SD{2n+l) ^ j dpC = -Mn J dpp^ I + y' {pf - {2T,a'f A[{pf . (3.4) 

Here, Nn = -^/^z)(2n+i)^S'" ^'^^ ^5" is the volume of the n-sphere. Since the action (|3.4|) 
does not depend on y{p) and At{p) explicitly we find two conserved charges 

1 i. (3.5) 



Sy'(p) (2,ra')SA[(p) 

General solutions at zero temperature Equations ( ^.5| ) can be solved in terms of y'{p) 
and A[{p) 

y'ip) 



1 d 



which can be integrated to (written in terms of hypergeometric functions 2^1) 

Note that we have chosen the integration constants such that y{p = 0) = and At{p = 0) = 0. 
The embedding functions behave near the AdS-origin (the zero temperature horizon) p — )• 
as 

yip) = -^^==p + 0{p'-^'), (3.9) 

for ah D3/D(2n+l)-systems. 

Boundary asymptotics From the asymptotic values limp_..oo y{p) = M and limp_!.oo ^tip) = 
p/ {2Tra') we get the relation between the constants of motion c, d and the physical parameters 
M and p foi n = 3: 

c = C3M {p^ - M^) , d = C3P (/x^ - M^) , (3.10) 

for n = 2 

c = - M2 , d = C2pVp'^ - ■ (3-11) 

where we have defined 

^ ~ r(i/3)3r(7/6)3 ' ^ ~ r(i/4)2r(5/4)2 ' ^^'^^^ 
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On-shell actions In order to extract thermodynamic quantities, we need the on-sheh ac- 
tion. Since the action ( |3.4D is divergent for p — )• oo, we introduce a cutoff A which we send 
to A —7- oo at the end of the calculation. To obtain a finite action we add appropriate 
counterterms of the form 

^counter = Ma dpp^ . (3.13) 

Thus the renormalized on-shell action density for the D3/D7- and D3/D5-systems is given by 

STen,n — •SD(2n+l) "l~ ■^counter • 

(3.14) 

4. Hyperscaling violation from D3/D7 and D3/D5 

In this section we briefly review the D3/D(2n-|-l)-brane setups at zero temperature, and the 
implications of their stationary analytic solutions for the thermodynamics of the dual theory 
[p^ ]. This will lead to the relations between the non-relativistic energy density e and the 
pressure p. 

4.1 Thermodynamics 

D3/D7 Since the density of the thermodynamical potential Vt in the grand canonical en- 
semble is given by the renormalized action density Sr-en; we obtain 

O = -Sren = _ liiZ^llW) (^2 _ ^2^2/3 ^ _]_c,{p^ _ ^2)2 _ (4 ^) 

4 



Note that the density of the thermodynamical potential is (minus) the pressure p, 

P=^C3(/x2-M2)2. (4.2) 

In order to obtain the density of the free energy, /, we Legendre transform fi, 

/ = n + /id= ^C3(/i^-M2)(3^2^M2). (4.3) 

Note that the density of the free energy, / is related to the (relativistic) energy density e by 
/ = e — sT, where s is the entropy density. Since we are working at zero temperature / = e. 

Let us now consider a special case, where /x is only slightly larger than M. Therefore we 
write p = M + fi with p <^ M. In this case we get to leading order 

n = -p = -CsM^fi^ , 
e = f = 2C3M3/2, 

c = d = 2C3M2/2 . (4.4) 

Note that in particular e ~ dM to leading order in p/M. In the following we consider p, as a 
non-relativistic chemical potential. The non-relativistic energy density e is then given by the 
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difference of the relativistic energy e and dM, i.e. e = e — dM where dM should be viewed 
as a rest-frame energy density, 

e = e-dM = C3M^fl^ , (4.5) 
to leading order in Jl/M. Note that to leading order in Ji/M we get 

e = p, (4.6) 
which concludes the proof of our first result in section |2[ 
D3/D5 Similarly we obtain 

^ = -Sren = - ^^^i^^^^\ d' - C^f' = -\c^{i^^ - M^f'^ . (4.7) 

The free energy density is then given by 

f = \C2^J^X''-M^{2^i^ + M''). (4.8) 
Considering also the limit fi = M -\- p, where p, /M ^ 1 we obtain to leading order in fi : 

e = f = V2C2M5/2/ii/2 ^ 

c = d = \/2C2M3/2^i/2 (49) 

Note that as in the case of D3/D7, also e ~ dM to leading order in fi. The non-relativistic 
energy e is given by 

/2 

e = €-dM = —C2M^/^fL^/^ , (4.10) 
3 

to leading order in p,. Comparing the non-relativistic energy density to the pressure we find 

e = f, (4.11) 

to leading order in p/M. 
4.2 Normal sound 

Here we compute the speed of sound. Its scaling with p will allow us to compute the values 
of exponents 9 and z. The speed of sound in a relativistic system is defined as 

vs' = '^. (4.12) 
oe 

Prom equations ( |4.3| ) and ( |4.4| ) we know that for the D3/D7-system we have p oc (^u^ — 
M2)2 and e tx (/U^ - M^)(3/i2 + M^). Therefore we can write 

, ^ {dp/df^ ^ jj^^-M^ ^ p_ 

(de/dfi) 3/x2 - M2 ^ M ' ^ ' 
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where in the last step we approximate p, <C M. For a propagating sound mode with frequency 
io and momentum k the scahng of the speed of sound in ( [4.13| ) imphes w = kyJji/M+0{k'^,p?) 
at small p, <^ M, and small uj <^ fl, k <^ p. This concludes the proof of our second result in 
section |2[ 

Let us repeat this computation for the D3/D5 system. We obtain 

where in the last step we again approximate p, <^ M. 

In fact, one gets Vg ~ p/M (and hence z = 2) in any scale invariant non-relativistic field 
theory (with or without hyperscaling violation) at zero temperature. As we have seen, in 
such a theory the free energy (and hence the pressure) scales as p^ with x = {n + z — 6)/z. 
The value of x will not be important for our argument. In a non-relativistic theory the speed 
of sound is given by 

2 'i- dp dp 

where p is the mass density, p = dM. The relativistic formula reduces to the non-relativistic 
one when we assume that e w dM, that is the energy density is dominated by the rest mass. 
Due to the scaling form of p we can then calculate as above 

, ^ J_ (dp/dp) ^ p_ 

' M{dd/dp) M' ^ ' 

for any x where a; 7^ 0, 1. 
4.3 Zero sound 

Now we make use of the dynamics of fluctuations on the D7-brane. In particular, we find 
the dispersion relations of the sound mode in these fluctuations, which was identified as zero 
sound in [^]. In our zero temperature, finite density system there is only one sound mode, 
which can be equivalently viewed as zero sound or normal sound. We will confirm that the 
dispersion relation of the probe brane sound mode is in perfect agreement with the normal 
sound in the previous subsection. For the standard formulas for the speed of normal sound in 
terms of thermodynamic quantities to hold we basically have to assume that the low energy 
dynamics of the system is governed by hydrodynamics. The reason zero sounds is usually 
viewed as different from normal sound is that the zero T, finite p regime traditionally is 
not identified with a hydrodynamic regime. Experience from holography however shows that 
for a hydrodynamic description having p to set the scale of what constitutes low energies is 
completely sufficient. We find it reassuring that the explicit fluctuation calculation agrees 
with the prediction from thermodynamics. Here we present a very condensed version of the 
computation. Details of the zero sound calculation can be found in appendix 

In order to calculate the zero sound, we consider fluctuations r/ of the form y(t, x, p) = 
yo{p) + v{'t, X, p) around the D7-brane embedding yo{p)- Furthermore, we choose the following 
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ansatz for the gauge field 

At{t, X, p) = At{p) + at{t, X, p) , (4.17) 

Ax{t,x,p) = axit,x,p) , (4.18) 

Ap{t,x,p) = ap{t,x,p) , (4.19) 

Ai{t,x,p) = 0. (4.20) 

Let us denote 

ftx{t,x,p) = dtax{t,x,p) - dxat{t,x,p) , (4.21) 

ftpit, X, p) = dtap{t, x, p) - dpat{t, x, p) , (4.22) 

fxpit, x,p = dxap{t, x, p) - dpQxit, X, p) . (4.23) 

From now on, we set R = 1 and 27ra' = 1. Then we can write the action to second order 
in the fluctuations a{t,x,p), r]{t,x,p) in the notation of [23] 

SdU+i) = -y / / '^Paip) [flp - f^(p)ftp - h{p)fl - hipMvf 

+h{p){dxvf + h{p){dpV? - 2h{p)ft.dxV 
-2Mp)ftpdpri] . (4.24) 

In order to generahze these expressions later, we have used the functions g{p) and fi{p)- For 
the D3/D7 brane model these are given by ( [AID with n = 6. 

After a Fourier transformation to frequencies w, momentum k (see (|A.2| )), change of 
variables from p to z = 1/p, and the field redefinitions X = krj and E = kat + u}ax, we obtain 
the equations of motion 

We solve these equations in two distinct limits. For the first case we take the limit z ^ 1 
of ( [4.25 ), find the solution to those approximate equations, and then expand that solution for 
ujz ^ 1. In the second case we take first the limit ujz <C 1, find the solution, and expand this 
solution for z 3> 1. Then we compare the low frequency, large z expansions of the solutions 
in both cases. From this matching we obtain the dispersion relation as 



which is easily expanded as w ~ k\JiilM at small <C M. This once again concludes a proof 
of our second postulate in section [2[. We explicitly determine the subleading coefficient /Ss 
in ( [OT[) . 
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Carrying out the zero sound computation for the D3/D5 system (see appendix we 
obtain the zero sound dispersion 




k + /32fc2 + . . . , (4.27) 



which again can be expanded as w k^/2fi/M at small /x ^ M. This confirms the scaling 
found from our computation of the normal sound speed ( 4.14| ). Again, we explicitly determine 



the subleading coefficient /32 in (|A.42[) . 

This completes the argument made in section |2|. The sound dispersion relations fix z = 2 
and then thermodynamics for both n = 3 and n = 2 gives 

z = 2, e = l, (4.28) 

which shows that hyperscaling is violated in the D3/D7- and D3/D5-probe brane systems 
near their phase transition. This is the main result of the present paper and we return to 
discussing it in section ^. 

5. Generalization to finite temperature 

In this section we turn on a small temperature T <^ ft. In order to derive the free energy we 
have to embed the massive D(2n+l)-brane into the background of Nc black 3-branes. Such 
massive D(2n+l)-brane embeddings at non-zero temperature and chemical potential are not 
known analytically. The blackening factors in the metric lead to more complicated differential 
equations for these embeddings compared to the previous case. 

Fortunately, for small temperatures T <C /i we can bypass the whole calculation. Ac- 



cording to 1 24] the correction of the free energy density to first order in the temperature T is 
given by 

/(/i, M, T) = fifi, M, T = 0) + 7TdT + 0{T^) . (5.1) 

where f{fi,T = 0) is the free energy of the relativistic system at zero temperature given by 
equations ( |4.3D and (^^) in the case of D7-branes and D5-branes, respectively. In order to 
obtain the non-relativistic result for the free energy we have to subtract the mass density dM, 
i.e. 

fnon-rel.i^^, M, T) = /(^, M, T) - dM . (5.2) 

In the case of D7-branes, the free energy to first order in the temperature is given by 

fnon-red^i, M, T) = hz^ii^^ " M^){^i? + + At^ ^iT - AM fl) + 0{T^) , (5.3) 
and in particular for fi = M + p, to leading order in fl/M, 

fnon-relXfi, M, T) = CsM^ fl^ 1 + ^ . (5.4) 

L /i J 
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We see that fnon-rei.{t^-,T) can be written in the form 

fnon-rel.ifi,T) fl^g , (5.5) 

where ^(0) = 1. In the analysis above we have determined the first coefficient in a Taylor 
expansion of the function g. Comparing with (O) we conclude that 

ziy = l, a = 0. (5.6) 

Since z = 2 we conclude z/ = 1/2. Together with = 1 and n = 3, we see that the modified 
hyperscaling relation 

{n + z-e)u = 2-a. (5.7) 

Let us repeat the analysis for the D3/D5 system. Then the free energy density to first 
order in the temperature reads 

fnon-redn.T) = ^C2M^'^JX^I^ + ^""^ + ^ (^)) " ^^'^^ 

In this case we see that 

zv = l, a = 1/2. (5.9) 

Note also that the assignment a = 1/2 = z = 2,n = 2 and 6 = 1 satisfies the modified 
hyperscaling relation. 



6. Discussion 

In this paper we have shown that the neighborhood of the second order phase transition 
governing the onset of finite density in the D3/D(2n+l)-systems is described by a non- 
relativistic, scale invariant field theory, which exhibits hyperscaling violation. The fixed point 
theory governs the non-relativistic limit where the chemical potential is just barely above the 
mass M. The quantities that scale non-trivially are the non-relativistic energy density (that 
is the relativistic energy density minus the rest mass of all the particles) as well as the non- 
relativistic chemical potential p, = — M . It is the appearance of this extra scale M that 
allows the free energy to scale differently than what one would expect from naive dimensional 
analysis. In both holographic systems we studied, one with n = 3 and one with n = 2 spatial 
dimensions, we find 9 = 1 and z = 2. To put our results into context, we should compare 
and contrast our strong coupling results with the same systems at weak coupling. At weak 
coupling, one would expect a finite chemical potential to drive the scalars into a condensed 
phase. The weakly coupled system should be well described by mean field theory. That is, we 
simply need to study the value of the effective potential at its minimum. Writing, following 
Landau, the most general potential for the scalar condensate |Xp consistent with symmetries 
one finds (in any spatial dimension n): 

V = {M'^ - /i^)|Xp + g\X\^ + . . . ~ -flX^ + gX^ + ... . (6.1) 
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Solving for the condensate one finds |Xp ~ (^u^ — M^) ~ fj, and hence the density of the 
grand canonical potential at the minimum scales as 

p ~ - M^f ~ fL^ . (6.2) 

We have shown that an equation of state of the form p ^ fi^ always gives rise to a sound 

_ n + z-9 

mode with dynamical critical exponent z = 2. Since we argued that p scales as p ^ ^ 
we see that p ^ p,'^ then fixes the hyperscaling violating exponent in mean field theory to 
6 = n — z. This mean field scaling is exactly what we find in the D3/D7 system with n = 3. 
For the D3/D5 system however mean field theory would predict 9 = while the brane system 
realizes 9 = 1. We find that a generalization of our scaling limit to small temperatures yields 
critical exponents a = and = 1/2 for the field theory dual to our D3/D7-system, as well as 
a = 1/2 and = 1/2 for D3/D5. With these values the modified hyperscaling relation ( |1.2D 
is satisfied for both of these theories. 

Perhaps even more interesting is the fact that for the D3/D5 system at zero temperature 
we have 9 = n — 1 which, according to the analysis of Q is exactly the value at which the 
system exhibits a logarithmic enhancement of the entanglement entropy as first pointed out 



in [25|. This behavior was seen as evidence for the existence of a Fermi surface in this system. 



It would be very interesting to see what the entanglement entropy is for the D3/D5 system. 



Unfortunately in order to use the Ryu-Takayanagi prescription [26] one seems to require the 
fully backreacted solution, so such a calculation is not currently feasible. Curiously, if we were 
to allow a potential of the form V ~ —JlX^ + g&X^ , then the mean field theory prediction 
would be = n — 1, consistent with our D3/D5 result. 

Our D3/D5-brane system realizes a field theory in its critical dimension, i.e. n = 2. 
Recently, a lattice computation of (/)^-theory at finite chemical potential in n = 3 dimensions 
has been performed [^]. We suggest to perform a similar lattice computation of the system in 
its critical dimension, namely 2 spatial dimensions, n = 2. It would then be very interesting 
to qualitatively compare our critical scaling results for the D3/D5-brane system in its critical 
dimension to such a lattice computation of 0^-theory in its critical dimension. Apart from 
bosonic matter the field theory dual to our D3/D5-brane system also contains fermions, whose 
influence on the critical exponents is not clear at this point. 

In this paper we have not considered the third D3/D(2n+l)-brane system with 4 ND 
directions, namely D3/D3 ||2^. Complications arise in this case from the fact that the bound- 



ary expansions of fields contain logarithms, see for example P9| , |30| , |31| . We will consider this 
case in a separate work, but already collect our zero sound results for D3/D3 in the appendix. 

The non-relativistic fixed points whose scaling behavior we have uncovered should clearly 
be further investigated. One could for example switch on external magnetic fields. The 
presence of a second scale would allow for non-trivial tests of our proposal. We have only 
analyzed in detail the 3 simplest probe brane system based on an AdSs x background. 
It would be very interesting to see if the general Dp/Dq system exhibits a similar behavior. 
In that case already the strongly coupled gauge theory itself violates hyperscaling Q. So it 
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would be very interesting to see what, if any, scaling is preserved at the non-relativistic fixed 
point. 
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A. Zero sound of D3/D7, D3/D5, and D3/D3 

Here we explicitly carry out the computation of zero sound modes in the D3/D(2n+1) systems 
studied in the body of this paper. Note that we follow in large parts the methods presented 



in m, |2|. 



In general the functions g{p) and fi{p) for a D3/D(2?i+l)-system with n = 1, 2, 3 are 
given by 



2 



(p2 + 2/o(p)2)2' ^^^^^ (p2 + yo(p)2)2' 
f t ^ y'o{p)^'t{p) f ( . y'o{p)K{p) 



{p' + yo{pW l + y',{pY-A',{pY 

f^^P) = p2n +d^-c^ 

The general background solutions are given by (3.6). 



Now we show how to obtain the fluctuation equations of motion in the form ( [4.25 ). In or 



der to express the field equations in a convenient way, let us perform a Fourier transformation 
to momentum space 

(2^)"+i g'^"""«Af(fc^,p), (A.2) 
___,ik.^-^^k,^py (A.3) 
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Using rotational invariance of the n spatial directions, we can consider without loss of gener- 
ality = {—ui, k,0, . . . , 0). The equations of motion in the gauge = read 

dplgifidpat - fjdpv)] - w kgf2ax - k^gf2at + k'^gfeV = , 

dpigdpaa:] + ujkgf2at + u?gf2ax - ^^kgfer] = , (A.4) 



dp[gihdpV + hdpat)] + 



uj^grj 



We also have to impose the Gauss constraint given by 

ka'x + ^fia't - iofyr]' = . 
Using the gauge invariant combination 



k'^gf'i'n - kujgfeax - k^gf&at = . 



(A.5) 



E = kat + uax , 

as well as the Gauss constraint, we can express a[ and a'^ in terms of 

k ^ "2 



A;2 - oj^f 



E' + 



kuj 



k^-oj-^fi k'^-UJ^f: 

Using the Gauss constraint, the first and second line of 
with two differential equations which look like 



-f7V'- 



(A.6) 



agree. Therefore we are left 



Jl 



g fi A;2 - a;2/ 



-) E' + fsf3{oJ^fi-k^)E- 



(A.7) 



f^kk"^ — huj"^ 



g h A;2 - a;2/ 



g h A;2 - uj'^fi 



E' + hh{uj^h-k^)E 



+ 



h k^-uj^fi k'^-uj^h 



^r]' + fsf3ioJ^fi-k^)v = 0- 



Here, we have used 



2^2 



Let us introduce 



_ /? + /l/5 _P^" 

h 

X = kr], 



p2n _ 



(A.8) 
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and change variables from p to z, given hy z = 1/ p. The equations of motion read 



E+\ - + - + ^^+ ■'^ 



z g fi k'^ - UJ^fi 



fi 



z g h A;2 - tj2y^ 



E+ P I g I ^7 ^ ^'^1 
\z g fj A;2 _ ^2y^ 



+ 



/s/; 



(A.c 



By adding and subtracting the equations of motion (with adequate prefactors) we can rewrite 
the equations of motion in the form given in ( 4.251) . 

We solve the two equations ( |4.25| ) in two different limits. First we consider the near- 
horizon limit z ^ 1. Using the precise form of /j, we can write the equations for the cases 
n = 3, n = 2 and n = 1 (D3/D7, D3/D5 and D3/D3 systems respectively) in the form 



2 (f -c^ 

E{z) + -E{z)+uj^ ?^E{x) = 0, 

z 



X{z) + -^X{z)+u;''' 



d? 



X{x) = 0. 



(A.IO) 
(A.ll) 



Note that the equations of motion are not coupled in the near- horizon limit. Using 



d2 ' 

we can write the solution of these two decoupled equations in the form 



E{z) = A- 



„iQz 



X(z) = B- 



(A.12) 



In addition, if we consider z in the limit Vtz <^ 1 



E{z) = — + iQ.A + . . . 



X{z) = — +inB + . 
z 



(A.13) 



Now consider ( [1.25 ) for all z but in the low frequency limit. In this limit we can drop 
the terms proportional to E and X. Therefore the equations read 



\z g hf^k'^-uj'^fi 
X +{-+-+ 



E 



hh k^ - uj^h 



■X = 0, 



z g hf^ A;2 - uj'^fi 



hh k^ - '^Vi 
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Let us now solve these equations of motion. It is convenient to write E{z) and X{z) in 
the form 



E{z): 

X{z) 
Then we obtain 

E{z) ^ 



1 nz) 

J. G{z) _ z^^-^Gjz) 

^ g{z) 22(l + (d2-c2)22n) ~ (1 + ((^2 _ c2)22n)3/2 " 



Z2(l + (d2_c2)^2n) 



_^ _ 2 _ 2n(cZ2 - c2)22n-l 

~5 ~ i ^ ~ 1 + (d2 - c2)z2n) 



(A.14) 

(A.15) 
(A.16) 



as weh as a similar expression for X{z) where we only have to exchange F by G. The equations 
of motion read 



F + 



hh - wVi 1 + (^2 - c2)z 



2n 



^ , / /5 fc^-CU^ 2n(d2-c2>2n-l \ ^ 



hf^k^-oj^fi l + {cP-c'^)z 



2n 



hh k^ - co^fi 

fi k' 
hh k"" - co^fi 



G{z) = 0. (A.18) 



After some rewriting, we see that these equations can be solved by the Ansatz 

F{z) = Di + 1)2^'" , 

G{Z) = Ci + C2^2n ^ 

we can determine D\ and D2 as a function of Ci,C2,uJ,k and c, d. We obtain 



Di = 
D2 = 



c2fc2C2 + (c2 - d'^){Cid'^ - G2)UJ^ 

cdk"^ 



Note that these two equation hold for all three cases n = 3, 2, 1. Then we have 

1 



F{z) 



cdk'^ 



C2 {k^ (-1 + C^Z^") + (1 + (-C^ + d^) Z^^) iJ") 

+Ci {J'k^ - (1 + (-c^ + d^) z^") u?) 



dxx^-^F{x) 



Therefore E{z) and X{z) are given by 

= Co + ^ (1 + (d2 _ c2)x2«)3/2 ' 

X{z) = Do+ ^ ^ 

Jo 



(1 + (d2 - c2)x2«)3/2 



(A.19) 

(A.20) 
(A.21) 



(A.22) 



(A.23) 



(A.24) 
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Note that for z — )■ 0, E{z) and X(z) behave as 

= Co + O(z"-1), X{z) = Do + O{z''-^) forn = 3,2, (A.25) 

and m the case n = 1, we have 

E{z) = Co + Ci\ogz + ... , X{z) = Do + Dilogz + . . . forn = 1 . (A.26) 

Since we want to determine the dispersion relation of the lowest-lying mode, we can set 
the sources Co and Dq to zero, i.e. Co = -Do = for the cases n = 3, 2. In contrast 
to that for D3/D3 with n = 1 the sources are not identified with the constant terms but 
with the coefficients in front of the logarithm, see also ||2^, 3C, 31|. So for D3/D3 we have 
Di = = Ci. Comparing this to ( |A.22| ) we immediately learn that our mode has to have the 
dispersion relation 

oj = k forn = 1 , (A.27) 

implying D2/C2 = d/c. Note that this is in agreement with ( A.13| ) if A, B, Co and Dq are 
chosen accordingly. 

Now let us perform the integration in equation ( |A.24 ) - note that F and G are given by 
the constants Ci and C2. For large z, in the cases n = 3, 2 we obtain 

E = Co + hC, + hC2 + ^ + ^ + ..., (A.28) 

z z 

X = Do + hC,+hC2 + ^ + ^ + ..., (A.29) 

z z 



while for n = 1, we get 



E = Co-—{d^-c^)-'/^ + ... , (A.30) 



z 



X = Do-^{d^-c^)-'/^ + .... (A.31) 

We obtain all the oi, 02, 61, 62 and ai, a2, 61, &2 explicitly, but refrain from listing them here. 

Now let us set Co = Dq = and let us match these solutions to the near-horizon solution 
in the cases n = 3, 2, i.e. 

61C1 + b,C2 + ^ + ^ = - + ^nA, (A.32) 

z z z 

61 Ci + 62C2 + ^ + ^ = - + . (A.33) 



Therefore 



feiCi + 62C2 = iO^, (A.34) 

aiCi + a2C2 = A, (A.35) 

biCi + 62C2 = inB , (A.36) 

aiCi + 02C2 = B , (A.37) 
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which means 



hCi + 62C2 - in{aiCi + asCa) = , 
hCi + 62C2 - in{diCi + asCs) = . 



(A.38) 
(A.39) 



In order to have a non-trivial solution for Ci, C2, the determinant of the matrix has to vanish. 
Inserting now co = ak + f^nk'^ + ■ ■ ■ we can determine a and /3„. The latter Ansatz for 00 is 
justified if we assume an analytic dispersion relation, non-integer powers of k would violate 
this assumption. 

For D3/D(2n-|-1) we finally get the dispersion relations 



k ~\~ Pnk , 



M2 - n/i2 

where the subleading terms are given for n = 3 by 



P3 



and for n = 2 by 
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